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I Abstract We present a new model universe based on the junction of FRW 

■ to flat Lemaitre-Tolman-Bondi (LTB) solutions of Einstein equations along 

QQ I our past light cone, bringing structures within the FRW models. The model is 

, assumed globally to be homogeneous, i.e. the cosmological principle is valid. 

C\J ' Local inhomogeneities within the past light cone are modeled as a flat LTB, 

, whereas those outside the light cone are assumed to be smoothed out and 

represented by a FRW model. The model is singularity free, always FRW far 
from the observer along the past light cone, gives way to a different luminos- 
ity distance relation as for the CDM/FRW models, a negative deceleration 
Qh' parameter near the observer, and correct linear and non-linear density con- 

trast. As a whole, the model behaves like a FRW model on the past light 
^ ' cone with a special behavior of the scale factor, Hubble and deceleration 

^ I parameter, mimicking dark energy. 



o 



o 



X 



Keywords Inhomogeneous Models • Dark Energy 



5— ( ' Paper in honor of Bahram Mashhoon's 60th birthday. 



Sh. Khosravi 

Physics Dept., Faculty of Science, Tarbiat Mo'alem University, Tehran, Iran. 
Institute for Studies in Theoretical Physics and Mathematics (IPM), Tehran, Iran. 
E-mail: khosravi@ipm.ir 

E. Kourkchi 

Department of Physics, Sharif University of Technology, Tehran, Iran. 
R. Mansouri 

Department of Physics, Sharif University of Technology, Tehran, Iran. 

Institute for Studies in Theoretical Physics and Mathematics (IPM), Tehran, Iran. 

Y. Akrami 

Department of Physics, Stockholm University, AlbaNova University Center, Stock- 
holm, Sweden. 



2 



1 Introduction 

Already in 1982, long before the detection of acceleration of the universe and 
suggestions in favor of inhomogeneous models as explanation for this novel 
cosmic effect, Bahram published a work in collaboration with Hossein Partovi 
[20] in which they developed a model with radial inhomogeneities as a first 
step toward a framework for comparison with observations. The paper did 
not get much attention due to the lack of observational evidence. The super- 
novae data indicating the acceleration of the universe [2311221 caused a revival 
of the idea of an inhomogeneous universe. This time the authors mainly con- 
sidered Lemaitre-Tolman-Bondi (LTB) solution of the Einstein equation as 
the model universe [T^[?7l l^ and indicated a change in the luminosity dis- 
tance relation towards explaining the acceleration of the universe due to the 
supernovae data (see [7] and [T6^ for references). Now, most of the authors 
consider the inhomogeneous solutions of the Einstein equations as an alter- 
native to FRW model of the universe (see [2],[H] and the references there). 
There is, however, another way of understanding the inhomogeneity and the 
interpretation of data on our past light cone. We know already from many 
simple examples that cutting a manifold along a null hypersurface and past- 
ing it again after some warping we may get gravitational shock or impulsive 
waves along the null hypersurface even in an otherwise Minkowski spacetime. 
This has been nicely shown in the paper by Penrose [21] . The cut-and-paste 
method is an elegant way to understand many astrophysical phenomena us- 
ing exact C° solutions of the Einstein equations. These solutions are usually 
everywhere except across a null hypersurface of the junction. Neverthe- 
less, the inhomogeneous models of the universe as partial explanation of dark 
energy, as Germans would say, are just beginning to be salonfdhig. The task 
force for dark energy [T] does not consider inhomogeneities as a possible 
explanation for dark energy. The first review or policy document in which 
the inhomogeneity is explicitly mentioned as a possible way to understand 
the acceleration of the universe seems to be "The Newsletter of the Topical 
Group in Gravitation ..." [TT] . 

We will report here on a model universe based on a cut-and-paste technol- 
ogy of matching an inhomogeneous model to a FRW model universe across 
the past light cone of the observer. The philosophy behind it is the following. 
After the decoupling era the universe, being almost in a global FRW state, 
begin to develop as a matter dominated universe witnessing a growth of struc- 
tures. The structures do not change the global homogeneous character of the 
universe, but locally the inhomogeneities begin to grow. Any local observer, 
in a globally Copernican view, will feel the effects of inhomogeneities in the 
course of time. Needless to say that the inhomogeneities outside his past light 
cone has no such effect and he may always assume that everything outside 
his past light cone is effectively homogeneous and is to be considered as a 
FRW model. On the contrary, the inhomogeneities in his vicinity within the 
light cone, which grows in time, may have direct cosmological consequences 
on his observations. These effects have to be increasing with decrease of the 
redshift on the past light cone. Therefore, we are faced with a model uni- 
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verse in which the past Ught cone of the observer is a boundary between a 
FRW model universe and an inhomogeneous model which goes to a FRW 
with increasing z, to match the requirement of almost FRW everywhere be- 
fore the onset of nonlinearities in the structures. For simplicity, we assume a 
LTB metric for the inhomogeneous part. The matching should be such that 
no shear is produced on the light cone as required for a simple boundary of 
observation. We will see that as a matching condition the light cone behaves 
as if we are in a FRW homogeneous model, but somehow dilated, mimicking 
a dark energy not existent in either FRW or LTB solutions on either side 
of the light cone. We have called this model Structured FRW (SFRW), as 
opposed to the familiar homogeneous FRW model of the universe [16]. 
The model we are considering is related to the Swiss Cheese model but differs 
from it substantially. The Swiss cheese model of cosmology, first suggested 
by Kantowski in 1969 [12], studies the effects of local inhomogeneities in 
the FRW models on the propagation of light through an otherwise homo- 
geneous and isotropic universe. Kantowski's model is constructed by taking 
the Friedman model {p = A = 0), randomly removing co- moving spheres 
from the dust, and placing Schwarzschild masses at the center of the holes. 
The remaining Friedman dust is interpreted as dust between clumps, and the 
point masses are interpreted as inhomogeneities. 

The paper has essentially two parts. In part one, we study in detail the flat 
LTB solution of Einstein equations, its singularities, and the bang time. This 
is just to get more insight how we are going to use the LTB which differs sub- 
stantially from its usage in cosmology in the recent literature. This is done 
in section 2 where we have written down the necessary formulas to define 
the LTB metric and their corresponding Einstein equations, the bang time 
and singularities, appropriate choice of the bang time, the past light cone, 
and the place of singularities within it. Part one may be skipped by those 
familiar with the LTB solution of the Einstein equations, except for the nu- 
merical results showing the place of singularities. In part two we define the 
model and study its cosmological consequences. Section 3 is devoted to the 
explicit definition of our SFRW model universe. In section 4 we look at the 
cosmological consequences of our model, the age of the universe, luminos- 
ity distance, deceleration parameter, and the density contrast, followed by a 
section on conclusions. 

2 Flat LTB solution and related cosmological quantities 

We constrain ourselves to the so-called flat or marginally bound LTB models. 
These are solutions of Einstein equations described by the metric 



in which overdot and prime (will thereafter) denote partial differentiation 
with respect to t and r, respectively. The corresponding Einstein equations 
turn out to be 



ds^ = -^c^dt^ + R'^ dr^ + R^{r, t){de^ + sin^ fld^^) 



(1) 



R^{r,t) 



2GM{r) 
R 



(2) 
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47rp(r, t) 



M'{r) 



(3) 



The density p(r, t) is in general an arbitrary function of r and t, and the 
integration time- independent function M{r) is defined as 



fRir.t) 4 

M(r) = 47r / p(r, t)R'^dR = — p(r, 



(4) 



where p, as a function of r and t, is the average density up to the radius R{r, t). 
The metric ([T]) can also be written in a form similar to the Robertson- Walker 
metric. The following definition 



/ , R(t,r) 
a{t,r) = 



brings the metric into the form 
ds^ = -c^df + a 



a r 
a 



dr^ 



(5) 



(6) 



For a homogeneous universe a doesn't depend on r and we get the familiar 
Robertson- Walker metric. The corresponding field equations can be written 
in the following familiar form 



^d-.2 _ G pc{r) 



(7) 



where we have introduced Pc(^) 



6J\/(r) 



indicating a quasi comoving r- 



dependent density. These are very similar to the familiar Friedman equations, 
except for the r-dependence of different quantities. The solutions to the field 
equations can be written in the form 



R{r,t) = 
a(r,t) = 



9GM(r) 



{t^tn{r))i. 



Now we choose the coordinate r such that 

47r 



(8) 



(9) 



where pc is a constant . To adapt this metric to observational data we need 
to know the backward light cone, the luminosity distance, the corresponding 
Hubble parameter, the deceleration parameter, the jerk, and the equation 
of state parameter w. We start with the radial light rays. The null geodesic 
corresponding to radially inward rays is given by 



cdt = —R'{r, t)dr. 



(10) 
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Assuming the redshift z as the parameter along the past hght cone, we obtain 

' (11) 



dz (l + z)7?'(r,i(r))' 



and 



dt 



-R'{r,t{r)) 



dz {1 + z)R'{r,t{r))' 



(12) 



where t(r) is evaluated along the rays moving radially inward according to 
pop. The luminosity distance is then given by 



Dl{z) = [l + zfR. 



(13) 



Definition of the Hubble function is not without ambiguity in the LTB mod- 
els. The reason is r and t dependence of the "scale factor a" [T71IT5] . Depend- 
ing on the use of the metric coefficients R or R' we may define Hj^ = ^ or 

Hr' — the Hubble parameter, which are different in general. We may 

also define a Hubble parameter as the expansion rate along the light cone 
|28j . It gives us a quantity which is easy to compare with the observations. 
Once we have the LTB-luminosity distance from (|13p . we interpret it as a 
distance in an effective FRW universe. Assuming the relation 



we invert it to define 



DLiz) = (1 + z) 



H{z) 



dz 
W) 

d DLiz)' 



dz 1 + z 



(14) 



(15) 



Note that in our structured FRW model this definition corresponds to the 
Hubble parameter of the FRW background outside the past light cone. There- 
fore, in general, we may have three different definitions of the Hubble param- 
eter. It has been shown, however, that in the case of glued LTB to FRW along 
a null hypersurface the three definitions coincide |13j . A numerical compari- 
son is given in section 4.4 and Fig. 14 
The associated deceleration parameter is then defined as 



-1 



1 + z 



H{z) 



dH{z) 
dz ' 



and the effective state parameter 



1/2] 



2(1 -I- z) d 



In 



H{z) 



(l + z)3/2 



3 dz 

In addition, we may define the jerk as in a FRW universe [251126] : 



J 



a a 



2iJ2 



(16) 



(17) 



(18) 



6 



where prime means derivative with respect to the argument a, and substitute 
an effective scale factor to obtain the corresponding value for our LTB based 
model. In this way we have a simple way to compare any LTB model with 
an effective FRW universe and also with the observational data. We may 
consider jerk as an alternative to w for parametrization of the dark energy. In 
our case it is to be considered as another way of interpreting inhomogeneities 
in terms of an effective FRW cosmology. 

We will see in section 3 that the effective FRW cosmological parameters 
defined here are the exact FRW quantities for our SFRW model on the light 
cone and that different Hubble parameter definitions for the LTB metric 
coincide with each other and are the same as the Hubble parameter of the 
background FRW on the light cone. 



2.1 Bang time and singularities 

The singularities of LTB metric, which are more sophisticated than in the 
case of the Robertson- Walker metric, have been discussed extensively in the 
literature (see for example [28], [S], [H]). Vanishing of each of the metric 

functions and its derivatives R, R', R, R, R', 

R" may lead to different singularities. In a general LTB metric there is an- 
other singularity, the event horizon, related to zero of 1 -I- where E{r) is 
the energy function of the LTB metric, absent in our flat LTB case. Vanishing 
of R"{t, r — 0) leads to the so-called central weak singularity in LTB models 
studied in [28] . For this central singularity in a flat LTB model to be absent 
one must have [28] 

C(0) = 0. (19) 

In selecting a bang time we will take this into consideration and will avoid 
those violating this condition. The place of the other singularities, irrespec- 
tive of its cosmological significance, are summarized as follows: 

R',R",R, R' ^oo 
R = 

R' = 
i?' = 

R" = 

The sign of t^, as it may be seen from the above relations, plays an important 
role in the discussion of the singularities. Taking note of these singularities, 
we will construct our model such that none of the singularities is on the light 
cone and all of them are far from the region where the LTB metric is effective. 



2 , 

1 , 

f ~ f rf 



(20) 



t = 



6t' + 3rt'' 
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2.2 Choice of the bang time 

According to SFRW model of the universe, the cosmological principle is valid. 
The universe from the point of view of each observer looks the same: inho- 
mogeneous locally but homogeneous far from the observer and everywhere 
outside its past light cone. Of course, outside the light cone we have ev- 
erywhere the local inhomogeneity which may affect us in future, but it is 
assumed to be smoothed out to a FRW homogeneous universe, as it is al- 
ways assumed in the familiar FRW model of the universe. The difference is 
on the light cone and in the vicinity of the observational point, where we do 
not want to assume the smearing out of inhomogeneities and would like to 
model it using a flat LTB solution of the Einstein equations, being glued to 
the outside flat FRW metric [13] . 

To formalize this requirement, we fabricate a bang time so that for r greater 
than a fiducial distance, say nL of the order of magnitude 100 MPc, the met- 
ric goes over to FRW, i.e. t„ — > const for r ^ nL. Here we assume L = 100 
MPc and leave n to be determined by the observation. For numerical cal- 
culation, we therefore define a dimension- less co-moving distance r' = 
For the sake of simplicity we rename from now on r' as r. Our coordinate 
r is now dimensionless and scaled by the inhomogeneity scale nL. But note 
that we may use r for the comoving coordinate or the scaled version of it 
interchangeably according to the context it is used! Now, for the bang time 
to have the desired property, we may write it in the following general form: 

p{r) + l + q{-) 

where p and q are polynomials in their arguments having no constant and 
linear term. The time factor a is another constant of the model in addition 
to inhomogeneity parameter n. It is obvious that for r ^ 1 the bang time 
approaches a constant, in fact zero for q 0, and for r 3> 1 it approaches a 
constant, zero again for p ^ 0, meaning that for large r we have effectively 
FRW metric again. Note that large r will corresponds in our case to redshifts 
bigger than 1. To see the cosmological effects of polynomials like p{r) we 
restrict it to powers of up to r^. On the other hand, it is easily seen that 
large powers of ^ have not much effect and we can ignore them within the 
scope of this paper. We, therefore, ignore them altogether to concentrate on 
the most significant cases and effects. Now, consider just the following two 
cases: 

t„ = - — % . (23) 

These bang times are just two examples of non-singular cases we will con- 
centrate on. In fact, there are many choices leading to almost FRW for large 
z, as we have tested numerically. They are, however, typically mal-behaved 
as far as singularities, density contrast, or deceleration is concerned, as it is 

2 

the case with i„ = "^^2 suggested in ^28j. We will see that the bang time 
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suggested in this paper is free of any singularity on the hght cone and reflects 
a reasonable behavior of the luminosity distance and the deceleration. 
We have plotted the behavior of t„ and tni as a function of r in Fig. [2l Both 
bang time functions decay very rapidly to zero, which is equivalent to rapid 
decay of the LTB metric to FRW. Therefore, it reflects the desired feature of 
SFRW. The tni, however, leads to singular behavior at the light cone for z- 
vlaues larger than 600. That is why we prefer to work with i„. In both cases, 
however, as it is shown in the following sections, the difference between LTB 
and FRW is almost negligible for redshifts z > 1. 

2.3 Past light cone and the place of singularities in the flat LTB model 

There are many papers dealing with singularities of the LTB metric. It is 
claimed that LTB metrics have either a weak singularity at the origin or do 
not lead to an acceleration [lOj . Here we report on the singularities of t„ and 
tni as defined above, just as two examples of bang time being non-singular at 
the origin. The question of acceleration is dealt with in the next section. In 
Figs. [3] and |4] we have plotted our past light cone and different singularities 
within it for t„ and tni- The singular points of R" are also sketched in the 
mentioned figures. This curve intersects the past light cone on a point which 
corresponds to a local maximum of R' as can be seen from Figs. [S] and 
[71 and none of the metric invariants, including that of Kretschman, have a 
singularity at this event. Therefore, we notice that all singularities are within 
the light cone, well before the time t — a, and well outside the vicinity of the 
light cone in the region where we have effectively FRW again. 

3 Definition of the model: Structured FRW Universe 

We are used to the Swiss cheese model as a random distribution of over- and 
under-dense regions in a constant time slice of the space-time, a constant 
time picture that does not take into account the realities of the observations 
along the past light cone and the possible evolutionary effects. We intend to 
modify this picture taking into account that all our observations are along 
the past light cone, and these observations are not affected by the events 
outside it. 

Let us accept the cosmological principle according to which all observers have 
the same picture of the universe. As the structures are effectively influencing 
the universe after the last scattering surface, our goal is to model the matter 
dominated and pressure-less universe. The universe before the last scattering 
surface is radiation dominated and structure-less, i.e. it is represented by a 
FRW metric. After the last scattering and with the growth of structures, our 
model gradually deviates from the simple CDM/FRW universe in the fol- 
lowing way. We will use a combination of FRW and LTB metrics connected 
to each other along the past light cone in a speciflc manner. Contrary to 
the familiar concept of using LTB metric to represent over- or under-density 
bubbles in the time constant slices of the universe in many papers published 
in the last years, we use it in a light cone adapted way. We will see that due 
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Present time 



(z=0| 



t4 




LSS (z-1100) 



Fig. 1 Observers A and B have the same picture of the universe: affected by the 
local inhomogeneities in their vicinity ( above the z = 1 plane); Outside their 
past light cone and far from the observer points A and B, the universe seems 
homogeneous. 



to our model construction the effect of LTB inhomogeneities vanishes on the 
light cone and we encounter another effective FRW which differs from the 
CDM/FRW, although we start with a CDM/FRW in combination with a 
pressure- less LTB. 

The universe is assumed to be globally FRW, but structured locally. Hence, 
any observer sees the universe in the following way: outside his past light 
cone, which is not observable to him, we assume the universe to be FRW, i.e. 
local inhomogeneities outside the light cone are assumed to be smoothed out 
and the effective metric is FRW. Note that the only use of this assumption 
is the induced Hubble constant on the light cone. Now, to implement the 
effect of the local inhomogeneities in the spirit of SFRW model, we assume 
the metric on the past light cone and its vicinity to be LTB. to define the 
vicinity we note that the order of magnitude of inhomogeneities is given by 
nL. Now at any z— value, the vicinity of the light cone can be determined 
according to the size of horizon. In fact, we are only concerned with the light 
cone behavior of the metric in this paper. In a following paper we will define 
the vicinity more precisely to study the averaging process of the Einstein 
equations and answer the question of the backreaction of inhomogeneities, 
its pitfalls, and how to make sense of it [14]. At any point outside this region 
the observer sees on the average, a FRW homogeneous universe again. Al- 
though, FRW and LTB could in principle be any of the three cases of open, 
flat, or close, it has been shown in ^13J that the only meaningful matching of 
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these two spaces along a null boundary is a flat-flat case. In addition, due to 
the cosmological preferences we will assume both metrics to be flat (Fig. [T]). 
Note the philosophy behind this picture: The universe is everywhere locally 
inhomogeneous, but homogeneous at large. Any observer, according to the 
cosmological principle, has the same picture. Now, he makes his model ac- 
cording to these rules: universe outside his past light cone, i.e. regions of the 
universe that has not influenced his past, is FRW. Locally, i.e. not far along 
the past light cone or for small z- values which will turn out to be of the order 
of z = 1, the observer sees a flat LTB on and within it. Far on the light cone 
and within it the model tends to be a homogeneous FRW universe again. 
This is of course in contrast to the usual assumption in the recent literature 
that the LTB represents bubbles of inhomogeneities. 
Now, let us assume the bang time to be 

tn = ^-^-T: (24) 

where r is scaled to nL. We will alternatively use the comoving coordinate 
r in the scaled form or not, and the reader may simply see from the context 
which one is meant. The constant a has the dimension of time. We have, 
therefore, two model parameters n and a to be determined by observation. 
Fig. [5] shows tn as a function of the redshift z, using pT|) . We have plotted 
r and i? as a function of z in the same diagram to see the coordinate, or 
redshift for which our LTB metric is essentially FRW. Obviously, the value 
z « 1, corresponding to r « 10 and t « 2a, is the boundary of transition 
from LTB to FRW. The inhomogeneity scale r ~ 10, or roughly 1000 Mpc, 
corresponds to the physical length « 1 Mpc at the time of the last scattering 
for z « 1100. Although the transition point is at a relatively small z, we 
can follow the effect of inhomogeneity up to the last scattering surface at 
z w 1100. The bang time is zero almost everywhere, except in our vicinity. 
Therefore, if the inhomogeneity has any effect, it must show up at our time, 
in accordance with the cosmological coincidence. 
Now, for this bang time ([M)) we have 

tn\r=O = 0. (25) 

Therefore, we do not expect any weak singularity at the origin [28] . In fact, 
for the LTB domain with this bang time, we have no singularity at all, as it 
is shown in Fig.[3l Vanishing of R" at the light cone and its vicinity indicates 
a maximum of R' for t — const, which reflects a feature of expanding layers. 
No invariant of the metric has a singularity within the domain of our interest. 
The singularity t = 3f„ which appears in the Kretschman invariant is well 
outside the region of our interest as shown in Fig. [3l We have also plotted in 
Fig. O i?' as a function of r for some fixed values of time to have a better 
understanding of variation of the LTB metric function. Note the behavior of 
i?' near its maximum values where R" vanishes and compare it with Fig. [3l 
The behavior of the LTB scale factor R/r and the metric function R' , for 
the sake of comparison the FRW scale factor for the CDM case (a = 0) as a 
function of z, is depicted in Fig. [TOl The density as a function of redshift is 
also plotted in Fig. [TT] and Fig. [12] for different ranges of z. 
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4 Cosmological consequences of SFRW and observational data 

We are now in a position to look at different cosmological consequences of 
the model. Angular distance and the luminosity distance are given implicitly 
by the equations dH] - [13]). We have to integrate from z = to the last 
scattering surface. The definition of the last scattering surface is a delicate 
problem. Its familiar value z « 1100 is based on a FRW model. But our 
SFRW model tends to a FRW at early times. Therefore, to match the SFRW 
at early times to a radiation dominated FRW we have to have the same 
age for the universe at the last scattering surface. This fixes the z-value in 
SFRW, which may differ from the one in FRW. It turns out, that for our 
choice of the bang time the difference is sensitive to the model parameters. 
Note, however, that for many other choices of the bang time there is no 
meaningful matching to the radiation dominated universe irrespective of any 
choice of the parameters. For the range of parameters we are considering, i.e. 
0.1 < n < 3 and a « 1 ( in scale of lO^^s) , the difference is negligible and we 
can fix the last scattering surface at z w 1100. Before going into its detail, 
we have to fix another constant of the model, namely pc- It is determined 
by the value of the scale factor of the universe at the last scattering surface. 
Taking it the same as that of FRW at the age of « 10^'^ s, we arrive at 
Pc ~ Wr^^gjcrr? or 10~'^^ in geometrical units. 



4.1 Angular distance and the age problem 

We have fixed now all the constants of the model. The last scattering surface 
is defined now by z = 1100, corresponding to the cosmic time t w 10^'^s. 
Integrating the p2)) . we end up with Fig. [13] for the light cone. For compar- 
ison, we have also plotted the CDM/FRW light cone for the same redshift 
interval. The age of the universe at the last scattering surface is almost the 
same for both models. Therefore, the value of time for z = determines the 
age of the universe in SFRW. It turns out to be 4/3 of the age in FRW, i.e. 
t = 3.9a = 15 X 10^ years. Therefore, there is no age problem in the model: 
the age of the universe in our model is well above any estimation we have in 
astrophysics and cosmology [5] (see Fig. [13]) . The angular distance defined 
by t) as a function of z is plotted in Fig. [9] 



4.2 Luminosity distance 

The model is now completely determined. The last scattering surface, having 
the same age in SFRW and FRW, corresponds to z « 1100. But as we have 
seen in the last section, the age of the universe is 4/3 of that in FRW. Now, 
we can integrate the equations pTl - [13]) to obtain the luminosity distance. 
Fig. [T51 shows the luminosity distance as a function of z. 
For the sake of completeness and more insight in the models under consid- 
eration, we have listed the results of likelihood analysis for the luminosity 
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distance of type la supernovae (SNe la - GOLD samples [24 ) given differ- 
ent forms of the bang time i„ in the table [1] ( assuming nL ~ lOOAfpc and 
a = 10" s). 



p{r) 


(for 182 SNe la) 


r* + l 


255 


r^ + r^ + 1 


265 


+ 1 


226 


+ l + r~'^ 


234 


+ r + 1 + r~'^ + 


239 


+ 1 + r"^ +r"^ 


251 


r + 1 


209 


r + 1 + r^^ 


220 


r + 1 + r"^ + r"^ 


226 


O.OOOlr'' +0.001r^ + 1 


286 




1700 


ACDM {Qa = 0.7) 


393 



Table 1 test for different choices of bang time. Bang time is defined as tn = 
and we put n — l,a — lQ^'^s,to — 3.9 x lO^^s. 



Different bang times t„ having a wide range of free parameters (i.e. a 
and n) enables us to get best fit luminosity distances, and likelihoods as in 
the other models like LCDM. Of course, we need to compare the results of 
our model with other cosmological data to choose the most suitable bang 
time and the acceptable range of its free parameters. At this stage, we found 
out that a should be approximately lO^^'s for any bang time tn- Looking at 
the density contrast, we can estimate the inhomogeneity factor as well (see 
section 4.4). Table[2]shows the best parameters for different bang times using 
SNe la data (GOLD samples): 



Banj 


I time 


a (xlO^'s) 


n 


(for 182 SNe la) 




+ 1)-^ 


1.1 


2.2 


188 


1 4 


+ 1)-^ 


1.0 


2.2 


212 


a(r 


+ r^ + l)-^ 


1.0 


2.9 


198 



Table 2 Best parameter choice for different bang times. 



4.3 Effective scale factor, Hubble parameter, and deceleration parameter 

There are different ways to define an effective scale factor for an LTB metric. 
Both R/r and R' could serve as effective scale factors, and correspondingly 

one may define ■§ or -§7 as the Hubble parameter. Of course, the definition 
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q(xIO^'s) n 

< 0.5 ~ Mot valid 

1 n ~ 1 S{zx 1100) ~ lO"'^ 

1 2 < 71 < 3 Valid for both conditions 
1.5 1 < n < 4 Valid for both conditions 

2 3 < n < 4 Valid for both conditions 
3<Q<5 4<n<5 Valid for both conditions 



Table 3 Density contrast for tn and different parameter values. 



(fT5|l we used to define the effective Hubble parameter may also differ from the 
other two definitions, ft is interesting, however, to note that all three defini- 
tions coincide if the LTB metric is matched exactly to FRW metric along the 
light cone, as has been shown in [13j . The requirement of the exact matching 
leads to a definite bang time, which could be calculated numerically. In fact, 
what we have done is the inverse problem: we have been looking for a suitable 
bang time mimicking the exact matching. What we have achieved is a bang 
time which leads almost everywhere, for the redshift values 1100 < z < 0.03, 
to the equivalence of all three Hubble parameter, i.e. to an exact matching. 
The actual problem of requiring exact matching and finding out the bang 
time is subject of a future work. 

Using relations and we have plotted in Fig. [T3] the different Hub- 
ble parameters. For comparison, we have also plotted iJ— values for CDM 
and LCDM models in the homogeneous FRW universe models (a = 0) . Note 
the peculiar z— dependence of the Hubble parameters for z < 0.5. There we 
observe a decrease of the Hubble parameter up to a value of about 50 and 
increase again to the present value of about 70. This behavior may be par- 
tially due to the deviation of Hfi from H = Hri . 

The effective deceleration parameter, as defined in ([16]), is plotted in FigfTBl 
for i„i, used in our SFRW model. It shows the effective deceleration for dif- 
ferent n— values. For z 1 all models almost coincide and have q = 1/2. 
They, however, differ from each other for small z. Negative deceleration is 
almost typical for all of them. Fig [TSl shows the deceleration parameter for t„. 



4.4 Density contrast 

Let us now look at the density contrast within the LTB domain. At any time, 
the density contrast within the LTB domain defined above is given by 

6= . (26) 

Calculation of S is straightforward. We expect 5 to be of the order of mag- 
nitude 10~^ for z w 1100 and of the order of magnitude 1 at the present 
time. Using t„, there is a rather wide range of a and n satisfying the den- 
sity contrast conditions. Table [3] shows the density contrast for some of the 
parameter values. 
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The above results for tn, and those of the luminosity distance {a = 10^^ s), 
shows that the inhomogencity factor n should be chosen in the range 1 < 
n < 3. Note that the density contrast should not be greater than the desired 
order of magnitudes. In the opposite case, one can compensate the lack of 
contrast by adding some density perturbations as it is done in FRW models. 
In the case of tn for a = 1 the inhomogeneity factor n has to be < 1 to 
achieve the correct density contrast at the last scattering surface. 

5 Conclusion 

The Structured FRW model of the matter dominated universe we are propos- 
ing is a singularity free model which incorporates the concept of local inho- 
mogeneities along the past light cone and is globally FRW. Although locally 
we have used the LTB metric, the effective model is a FRW type;: along the 
light cone different Hubble parameters coincide, except for small values of z, 
as if we have just a FRW model. Therefore, we could consider this SFRW 
as a modified CDM/FRW with novel features. This is in contrast to familiar 
use of LTB metric to model local over- or under-dense bubbles. Not only the 
age of the universe is increased, the deceleration parameter has also an in- 
teresting behavior: it is negative for small z- values and goes to 1/2 for larger 
rcdshifts. The density contrast changes from the order of magnitude 10^^ at 
the time of the last scattering surface to 1 for the present time. Our model 
has two parameters a and n both of the order of magnitude 1. The bang 
time is defined in a way that there is no central weak singularity present in 
the model; although it is not unique, the difference between different bang 
times in the class we have defined is marginal. We conclude that allowing for 
local changes along the past light cone one may construct model universes 
that behave Friedmanian but differs from the CDM model and leads to local 
acceleration, at the same time increases the age of the universe and refiects 
the growth of density contrast. 
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ments. R.M. would also like to thank ISMO (Tehran) for partial support through 
TWAS-IC funds. 



References 

1. Albrecht, A., Bernstein, G., Cahn, R., Frcodman, W.L., Hewitt, J., Hu, W., 
Huth, J., Kamionkowski, M., Kolb, E.W., Knox, L., Mather, J.C., Staggs, S., 
Suntzeff, N.B.: Report of the Dark Energy Task Force, astro-ph/0609591 (2006) 

2. Biswas, T., Mansouri, R., Notari, A.: Nonlinear Structure Formation and "Ap- 
parent" Acceleration: an Investigation, astro-ph/0606703 (2006) 

3. Blanchard, A., Douspis, M., Rowan-Robinson, M., Sarkar, S.: An alternative 
to the cosmological "concordance model". Astrou Astropiiys 412, 35 (2003). 
DOI 10.1051/0004-6361:20031425 

4. Bondi, H.: Spiicrically symmetrical models in general relativity. Mon Not R 
Astro Soc 107, 410 (1947) 

5. Carretta, E., Gratton, R.G., Clementini, G., Fusi Pecci, F.: Distances, Ages, 
and Epoch of Formation of Globular Clusters. ApJ 533, 215 (2000). DOI 
10.1086/308629 



15 



6. Celerier, M.N.: Do we really see a cosmological constant in the supernovas 
data? Astron Astrophys 353, 63 (2000) 

7. Celerier, M.N.: Is the apparent acceleration of the Universe expansion driven 
by a dark energy-like component or by inhomogeneities? astro-ph/0612222 

(2006) 

8. Celerier, M.N.: The accelerated expansion of the Universe challenged by an 
effect of the inhomogeneities. A review, astro-ph/0702416 (2007) 

9. Christodoulou, D.: Violation of cosmic censorship in the gravitational collapse 
of a dust cloud. Communications in Mathematical Physics 93, 171 (1984) 

10. Flanagan, E.E.: Can suporliorizon perturbations drive the acceleration of the 
Universe? Pliys Rev D 71(10), 103,521 (2005). DOI 10.1103/PhysRevD.71. 
103521 

11. Garfinklo, D.: Matters of Gravity, The Newsletter of the Topical Group in 
Gravitation of the American Physical Society, Volume 29, Winter 2007. gr- 

qc/0702081 (2007) 

12. Kantowski, R.: Corrections in the Lumiuosity-Rcdshift Relations of the Homo- 
geneous Fricd-Mann Models. ApJ 155, 89 (1969) 

13. Khakshournia. S., Mansouri, R.; Matching LTB and FRW spacetimes through 
a null hypcrsurfacc. gr-qc/0702130 (2007) 

14. Khosravi, S., Kourkchi, E., Mansouri, R.: Modification of averaging process in 
GR: case study SFRW. In Preparation 

15. Lemaitre, G.A.: L' Univcrs en expansion. Ann. Soc. Sci. Bruxelles p. 51 (1933). 
Reprinted in GRG, 29 (1997), 641. 

16. Mansouri, R.: Structured FRW universe loads to acceleration: a non- 
perturbative approach, astro-ph/0512605 (2005) 

17. Moffat, J.W.: Cosmic microwave background, accelerating universe and inlio- 
mogeneous cosmology. Journal of Cosmology and Astro-Particle Physics 10, 
12 (2005). DOI 10.1088/1475-7516/2005/10/012 

18. Moffat, J.W.: Late-time inhomogeneity and acceleration without dark energy. 
Journal of Cosmology and Astro-Particle Physics 5, 1 (2006). DOI 10.1088/ 
1475-7516/2006/05/001 

19. Newman, R.P.A.C.: Strengths of naked singularities in Tolman-Bondi space- 
times. Classical and Quantum Gravity 3, 527 (1986) 

20. Partovi, M.H.. Mashhoon, B.: Toward verification of large-scale homogeneity 
in cosmology. ApJ 276, 4 (1984). DOI 10.1086/161588 

21. Penrose, R.: The Geometry of Impulsive Gravitational Waves. In: 
L. O'Raifcartaigh (ed.) General Relativity. Oxford Clarendon (1972) 

22. Perlmuttcr, S., Aldering, G., Goldhaber, G., Knop, R.A., Nugent, P., Castro, 
P.G., Deustua, S., Fabbro, S., Goobar, A., Groom, D.E., Hook, I.M., Kim, A.G., 
Kim, M.Y., Lee, J.C., Nunes, N.J., Pain, R., Pennypacker, C.R., Quimby, R., 
Lidman, C, Ellis, R.S., Irwin, M., McMahon, R.G., Ruiz-Lapuente, P., Wal- 
ton, N., Schaefer, B., Boyle, B.J., Filippenko, A.V., Matheson, T., Fruchter, 

A. S., Panagia, N., Newberg, H.J.M., Couch, W.J., The Supernova Cosmology 
Project: Measurements of Omega and Lambda from 42 High-Redshift Super- 
novae. ApJ 517, 565 (1999). DOI 10.1086/307221 

23. Riess, A.G., Filippenko, A.V., Challis, P., Clocchiatti, A., Diercks, A., Gar- 
navich, P.M., Gilliland, R.L., Hogan, C.J., Jha, S., Kirshnor, R.P., Loibundgut, 

B. , Phillips, M.M., Reiss, D., Schmidt, B.P., Schommer, R.A., Smith, R.C., 
Spyromilio, J., Stubbs, C, Suntzeff, N.B., Tonry, J.: Observational Evidence 
from Supernovae for an Accelerating Universe and a Cosmological Constant. 
Astron J 116, 1009 (1998). DOI 10.1086/300499 

24. Riess, A.G., Strolger, L.G., Casertano, S., Ferguson, H.C., Mobaslicr, B., Gold, 
B., Challis, P.J., Filippenko, A.V., Jha, S., Li, W., Tonry, J., Foley, R., Kirsh- 
ner, R.P., Dickinson, M., MacDonald, E., Eisenstein, D., Livio, M., Younger, 
J., Xu, C, Dahlen, T., Stern, D.: New Hubble Space Telescope Discoveries of 
Type la Supernovae at z > 1: Narrowing Constraints on the Early Behavior of 
Dark Energy, astro-ph/0611572 (2006) 

25. Riess, A.G., Strolger, L.G., Tonry, J., Casertano, S., Ferguson, H.C., Mobasher, 
B., Challis, P., Filippenko, A.V., Jha, S., Li, W., Chornock, R., Kirshncr, R.P., 
Leiburidgut, B., Dickinson, M., Livio, M., Giavalisco, M., Steidel, C.C., Bem'tez, 



16 



T., Tsvetanov, Z.: Type la Supernova Discoveries at z > 1 from the Hubble 
Space Telescope: Evidence for Past Deceleration and Constraints on Dark En- 
ergy Evolution. ApJ 607, 665 (2004). DOI 10.1086/383612 

26. Sahni, V., Saini, T.D., Starobinsky, A.A., Alam, U.: Statefinder-A New Geo- 
metrical Diagnostic of Dark Energy. Journal of Experimental and Theoretical 
Physics Letteres 77, 201 (2003) 

27. Tolman, R.C.: Effect of Inhomogencity on Cosmological Models. Proceedings 
of the National Academy of Science (USA) 20, 169 (1934) 

28. Vanderveld, R.A., Flanagan, E.E., Wasserman, I.: Mimicking dark energy with 
Lemaitre-Tolman-Bondi models: Weak central singularities and critical points. 
Phys Rev D 74(2), 023,506 (2006). DOI 10.1103/PhysRevD.74.023506 



17 



x10^'(s) 



- 




\ 


- 




\\ - 


- 


t - " 

" r4 + r2 + l 


t\ 
1 1 

\\ 


- 




\\ 










\\ 






l\ 
\\ 






l\ 
\ \ 






\ \ 
\ \ 






\ \ 
\ \ 






\ \ 
\ \ 



0' — — — — ■ ^ ■" ^ ' 

0.01 0.1 1 10 100 1000 10000 

r (MPc) 



Fig. 2 t„ vs r using different functions. The blue curve for t„ = ^i^"-2_^_-i and the 
red curve for tni = ^" . 
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Fig. 3 Light cone is shown (time vs io(?io(r)) for tni = r'"+i ' ('^ 
The sohd hne curve shows the sheU focussing area R = at t — t„. On this curve 
both R' and R" are singular. Here, one can easily see that one of these singular 
points lie on the light cone. This singular point maybe avoided if we restrict the 
value of n to be less than 0.2 or choose another bang time t„(r-). The shell focussing 
singularity crosses the light cone at z— values well above the last scattering surface 
which is outside the range of applicability of our model, although it may also 
be avoided using suitable inhomogeneity parameter. The dash-dotted green curve 
corresponds to Kretschman singularity (t = 3t„(r-)). 
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x10''(s) 



Singularity Diagram : t 
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Fig. 4 Similar to Fig. 3, for t„ — ^.i^'^^i^i , {ct ~ 10^^ s,n — 1). All singularities 
are well inside the curve t — 3t„(r). The R' singularity points (corresponding to 
t = tn — ^rt'„) are not shown because they almost lie completely on the curve for 
t = tn{r). 
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Fig. 5 B! is plotted versus r for tn and different constant times. Looking at i?' as 
an effective scale factor, it shows that the scale of the universe increases with time, 
although the rate of cosmic expansion is different at different places. 
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Fig. 6 R' vs r for t„ and t — 3 x 10^^ s. No singularity is seen here. 
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Fig. 7 A typical singular behavior of R' inside the light cone for t„ and t = 
8 X 10^'' s. The bump after the singular point corresponds to the vanishing of R" . 
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Fig. 8 A typical singular behavior of R' vs r for £„ and constant £ = 8 x 10^® s. 
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Fig. 9 tn, R, r is drawn versus z. The bang time is almost zero for z > 0.1 
. When we go back in time along the light cone, the scale of the universe (blue 
curve) goes to zero (i.e. ii — > 0). 

1.4 1 




z 



Fig. 10 The SFRW effective scale factor R/r as a function of z for t„ (blue curve). 
The red points show R' as a function of z. Both are almost the same for z > 0.1. 
The effect of inhomogeneities is more clear for z < 0.1 where R' 7^ R/r. The green 
curve shows the scale factor for a CDM/FRW model (a = 0). 
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Fig. 11 Density p as a function of z for tn in low redshift region 
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Fig. 12 Density p as a function of z for t„ in high redshift region 
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Fig. 13 Comparison of light cones for both SFRW (LTB) and FRW (CDM) 
models. The age of the Universe is greater in SFRW model by the factor 4/3. We 
got these results for z < 1100 (after last scattering time). 
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Hubble Parameter 




Fig. 14 Using the luminosity distance and ea llSI the effective Hubble parameter is 
plotted, ft shows also the Hubble parameter for CDM and ylCDM models. Compar- 
ison of the effective Hubble parameter based on the luminosity distance (Hsfrw) 
and R' /R' shows, for all values of redshift parameter, these two definitions are the 
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Fig. 15 fi^m-M vsz (redshift) for diflerent models (LCDM, CDM, t„ SFRW). 
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Fig. 16 The effective deceleration parameter for SFRW model {tni — 73^) using 
the effective H{z) and ea ll6l This diagram clearly shows that, for small values of 
2, q is negative and therefore, the universe is effectively accelerating at the present 
time. Increasing the parameter n, the transition of acceleration to deceleration 
occurs in higher redshifts. 
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Fig. 17 The effective state parameter for SFRW model {tni = -prpi) using the 
effective H{z) and egll?! 
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Fig. 18 Similar to Fig[T6] using t„ 
0.0001, B = 0.001. 
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Fig. 19 Similar to Fig[T3 using t„ — -p^ 



and t„ 



